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Understanding how new opinions spread through a community or how consensus emerges in noisy 
environments can have a significant impact in our comprehension of the social relations among 
individuals. In this work a new model for the dynamics of opinion formation is introduced. The 
model is based on a non-linear interaction between opinion vectors of agents plus a stochastic variable 
to account for the effect of noise in the way the agents communicate. The presented dynamics is 
able to generate rich dynamical patterns of interacting groups or clusters of agents with the same 
opinion without a leader or centralized control. Our results show that by increasing the intensity of 
noise, the system goes from consensus to a disordered state. Depending on the number of competing 
opinions and the details of the network of interactions, the system displays a first or a second order 
transition. We compare the behavior of different topologies of interactions: Id chains, annealed and 
complex networks. 
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I. INTRODUCTION 

An interesting application concerning the structure of 
social networks 0, Q is the modeling of the dynamics 
of opinion formation. Specific measurements that cha- 
racterize the statistics behind the existence of different 
groups and affiliations within human populations, justify 
to model such aspect of human behavior. The idea be- 
hind this field is to find simple rules of interactions behind 
the nodes or agents, each of which carries its own chang- 
ing color or opinion, trying to reproduce the emergence 
of complex patterns observed in reality. 

Such opinions can be defined by a finite number of 
integers as in the model proposed by Sznajd et al. [3| 
or can even be represented by real numbers, having a 
rich spectrum and opening the possibility for having as 
many opinions as agents; like in the model proposed by 
Deffuant et al. Q. In both cases the proposed dynamics 
has a natural absorbing state or consensus, in which all 
the agents share the same opinion. Other models as the 
one by Hegselmann and Krause @, the voter model 0, 
the Galam's majority rule Q , and the Axelrod's model [8| 
are reviewed in Ref. Q. As pointed out there, very few 
and non conclusive results exist for the consensus models 
on complex networks. 

In this work, we present a general model, where opi- 
nions are represented by vectors with real components 
and the agents interact with a non-linear rule. We pro- 
pose, for the abstract space of human opinions a dynam- 
ical rule where each agent has a opinion vector that is 



fixed in modulus. Every time step an agent interacts 
with its neighbors and assumes a new value for its opin- 
ion vector that is a function of the average direction of its 
neighboring agents plus an added noisy term. The result- 
ing behavior presents two important characteristics, (i) 
although the model allows a continuous change from one 
opinion to some other, the interaction favors extremely 
decided states over undecided states, and (ii) the system 
ubiquitously evolves to coordination and grouping with- 
out the need of leaders or centralized control. The fact 
that the modulus of the opinion vector is constant de- 
scribes the strength of an opinion about a specific topic 
at the expense of the other beliefs. According to our 
model, undecided agents (i.e., those that do not have 
a strong belief in one particular opinion) can not affect 
the ones with a strong opinion. This type of interac- 
tion is somewhat different to the one used in models of 
opinion formation, which usually consists in weighted av- 
erages. Similar rules have been proposed for example in 
Ref. [l0( to explain how very large populations are able to 
converge to the use of a particular word or grammatical 
construction without global coordination. 



We study the transition to consensus as a function of 
noise. As a new feature, we find that different types of 
transitions with or without hysteresis are observed de- 
pending on the dimension of the opinion vector. Addi- 
tionally, we observe that the transition is controlled by 
the interaction dynamics and is independent on the cor- 
relations of the network topology. 
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We start in Section |TT] by describing in detail the pro- 
posed model for interactions between opinions. In Sec- 
tion lllll we present the results of the transition to consen- 
sus as a function of the noise for Id chains, annealed, and 
complex networks. Conclusions are given in Section [IVI 
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FIG. 1: Finite size scaling of the transition in a system with 
a two-dimensional opinion vector with random interactions. 
The figure compares the frequency of agents with the same 
dominant opinion versus the noise r\ for different population 
sizes N. a: Original data, b: Finite size scaling. Each point 
corresponds to an average over 10 — 20 runs with different 
random seeds. 



II. THE OPINION MODEL 

The system comprises a fixed number of N agents. Ev- 
ery agent i is characterized by its own opinion vector 
<2n of n — 1, .., O opinions. Each element of this vector 
corresponds to a different opinion about the same topic. 
Negative values of the elements are not allowed. The sum 
of the vector elements of an agent is one: J2n=i = 1- 
For instance, agent j favors communism with 20% and 



capitalism with 80% (given O = 2): a{ 



0.2, a 2 W 



0.; 



Each time step every agent actualizes its opinion vector 
by comparing its values to the ones of its ki nearest neigh- 
bors. These other agents are chosen by the topology of 
the graph, and the agent updates its opinion vector due 
to the following rule, 



a«(i) = >>W(i)a«(i) + M*)> 



i=i 



(1) 



where g(t) presents a stochastic variable, distributed uni- 
formly in the interval [0, rf\ . With this exclusively positive 

noise we assure that a„ (t) > 0. This stochasticity can 



be interpreted to be due to misunderstandings among 
the agents, the spread of wrong information, or other 
perturbing actions. 

The interaction term in this model is of second order. 
Thus, in a noiseless environment, the agents tend to have 
the same stronger opinion. The factor ki avoids that 
agents which have more connections feel less noise. In 
order to guarantee that the sum of opinions is equal to 
one, the vector is normalized afterwards similar to the 
model presented in Ref. [Tl|, 
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In order to elucidate the principal properties of the up- 
date rule given by Eqs. (Q} and ([2]), we examine in detail 
the noiseless interactions between three types of agents 
with different characteristic values of a two-dimensional 
opinion vector (O — 2), namely, a\ = {0.8,0.2}, a 2 = 
{0.5,0.5} and as — {0.2,0.8}. First, an interaction be- 
tween an agent having ai with another having the same 
a\ results in {0.94,0.06} — interactions between agents 
with the same dominant opinion strengthens their belief 
in this opinion. a\ with a 2 yields {0.8,0.2} — interactions 
with "undecided" agents are ineffective in the sense that 
agents without dominant opinion are not able to convince 
another agent. On the other hand, this interaction will 
have a substantial effect on the undecided agent, i.e., 
undecided agents are convinced easily. The interaction 
between a\ and results in {0.5, 0.5} — interactions of 
agents with opposite opinions lead them to become less 
decided. 

At the beginning of a simulation the opinion vectors 
are initialized either randomly or by consensus: for ran- 
dom initialization we toss for each opinion component of 
each agent a number between zero and one. The opinion 
vectors are normalized afterwards according to Eq. 
The other way to initialize the system (consensus) is by 
setting one to the first element of each vector and fill the 
rest with zeroes. 

The main parameter of this model is given by the max- 
imal noise i] which we will call from now on the control 
parameter. Its role corresponds to the one of a tempera- 
ture in physical systems. In a social system, the noise 
represents any internal or external interference in the 
communication among the agents. Other free parame- 
ters of the system are given by the number of agents N, 
the number of opinions O and the number of agents ki to 
interact with per time step. The last parameter can be 
different for distinct agents depending on the topology of 
the actual network. 

A simple mean-field solution of the model without 
noise can be derived. Suppose a state where all agents 
have the same values in their opinion vectors. Thus the 
index of the agents can be suppressed, cin\t) = a n (t) 
and (t + 1) = a n (t + 1), and Eqs. (|T|) and @ can be 
summarized. In the case of two opinions the equations 
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FIG. 2: The transition of the system for different numbers 
of opinions O and different initial conditions. The popula- 
tion size is 500. a: The outcome of a system of two opinions 
(circles) is independent of its initialization, b: In the case of 
three opinions (dashed line), the curves present a hysteresis 
and the results are different, if the field is initialized randomly 
(squares) or with the system being already in the consensus 
state (stars), c: A system with an opinion vector containing 
ten opinions also exhibits hysteresis (dash-dotted line).d Av- 
erage value of the dominant opinion a 0mam vs. n for the same 
simulation as in c. 



are a sampling of the whole system, it is expected that 
this annealed approximation should behave similar to a 
me an- field. 
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FIG. 3: Ten opinions on a one-dimensional chain. These re- 
sults illustrate simulations without noise (jj — 0) and a system 
of 1,000 agents, a: The figure zooms on the first 200 agents 
of the population, where each line corresponds to a different 
element of the opinion vectors. Only the first five opinions are 
displayed. The agents form local clusters of different domi- 
nant opinions, b: The cluster sizes distribute following an 
exponential decay. 



correspond to the map, 
a\{t) 



ai(m) 



o?(t) + «*!(*) 



, a 2 (t+l) 



al{t) 



al(t)+al{t)' 



(3) 



The fixed points of these equations are (01,02) £ 
{(0, 1), (1,0), (0.5,0.5)} where the first two ones are sta- 
ble and the last one is unstable. The solutions for O opin- 
ions are in {(1, 0, 0, ...), (0, 1, 0, 0, ...),..., (0, 0, 0, 1)} 
with all a n stable. All other solutions have at least one 
unstable element of the opinion vector and thus the un- 
stable element influences the other ones until an absorb- 
ing state with one opinion totally dominant is reached. 



III. RESULTS 
A. Annealed interactions 

First, we present simulations of the model without 
fixed topology. Each time step, a simulation runs over 
all agents. For each of them and at each time step, two 
new random partners are chosen to interact. We chose 
interaction with two other agents (ki = k = 2) in order 
to facilitate the comparison of this case with the one of 
a one-dimensional chain which will be explained in the 
following section. The annealed approach avoids long 
term behavior and the distribution of opinions reaches 
the stationary state fast. Because the interacting units 



The results reveal that the system can reach two dif- 
ferent absorbing states. At small values of the control 
parameter (maximum noise rj), one opinion completely 
dominates the system, o max . For a noise rj larger than 
a certain value, each opinion remains with the same fre- 
quency, I/O. The order parameter D is the frequency of 
the agents which have an opinion vector with the same 
dominant opinion, being itself dominant in the system. 
More precisely: for each agent we search its strongest 
opinion and then count, for each opinion, the number 
"n" of agents with this opinion as their dominant one. 
The largest value no, and so the most dominant one of 
the system, determines D = n^/N. (D) means, that 
we average D over many time steps. This order param- 
eter is normalized, so that it is unity if all agents have 
the same dominant opinion, a state we call the consensus 
state. The value I/O corresponds to a uniform distribu- 
tion of opinions. A transition occurs between consensus 
and uniform distribution, when (D) goes from 1 to 1/2 
in the case of two opinions (Fig. QJi). 

The transition becomes more abrupt for larger pop- 
ulation sizes. A transition point characteristic for the 
jump from the consensus to the uniform states is located 
at n c w 0.5, increasing with the population size. This 
transition seems to be a phase transition of second order. 
We carried out finite size scaling in order to obtain the 
critical exponents (Fig. []J>). Near the critical point the 
curves coincide using the scaling relations, 
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FIG. 4: The size of the largest cluster increases in time until 
it reaches the population size TV = 1,000. Here, we see the 
results of different random initializations of the system on a 
chain. The number of opinions is ten, and the noise n = 0.2. 



with v w 2.4 ± 0.1, = 0.15 ± 0.05, and the critical noise 
r] c = 0.52. 

Fig. [2] shows that in the case of annealed interactions 
the transition becomes of first order for simulations with 
an opinion vector of more than two opinion elements. 
The fluctuations do not increase at the transition point. 
Now, the transition from the consensus to the uniform 
state depends on the initialization and is much more 
abrupt. If the initialization is random, the system jumps 
to the consensus state at lower values of 77 than in the 
case of an initialization with consensus in one opinion. A 
transition with a typical hysteresis occurs at lower values 
of r\ if we increase the dimension of the opinion vector. 

Note that D gives us the fraction of agents with dom- 
inant opinion (o max ) but does not contain information 
about a Gmax , the magnitude of the component associated 
to Omax- In Fig- [IJi we plot a 0maw vs. r/ for the same 
simulations presented in Fig.[2t. a Qmax is larger for lower 
values of 77 and below a certain r] c consensus is observed 
for both kinds of initializations, only when a large value 
of a ,„ ax is reached. This is a nice feature of our model: 
consensus and resolution emerge together in the system. 
That is, the agents can only make up their minds for a 
preferred opinion when consensus is achieved through all 
the system. 



B. One-dimensional topology 

If we put the agents on an one-dimensional lattice with 
periodic boundary conditions, or, in other words, a chain, 
the results become different. First, we concentrate on the 
case of ten opinions and no noise (77 = 0): The system 
is now highly dependent on the initial state. A random 
initialization of the opinion vectors leads to the situation 
depicted in Fig. [3^. The same amount of each opinion 



seems present in the system during the evolution. The 
system organizes itself by rearranging its opinion vectors 
to form local clusters of different sizes. In one cluster the 
same opinion dominates for all agents. Each agent has a 
well pronounced dominance of an opinion (its value being 
nearly one) , and the interfaces between clusters of differ- 
ent dominant opinions are very sharp. These clusters 
develop fast after the beginning of the simulation. The 
distribution of cluster sizes follows the exponential de- 
cay of a Poisson distribution (Fig. (SJa). The results with 
77 = are qualitatively the same for different numbers of 
opinions, O. 

Noise (77 > 0) leads to a slow increase of one of the 
ten opinions with time. The dominant opinion absorbs 
more and more of the losing opinions. Fig. [?] illustrates 
how the largest cluster of the system temporally evolves 
for 7/ = 0.2. As also can be recognized in this figure, 
the time to reach consensus can be really long, even in a 
small system of 1, 000 agents. 

With non-zero but small noise, the information prop- 
agates slowly through the sample. Because of that, the 
time to reach the absorbing state is much larger than 
in the case of random interactions. Next we consider 
a system consisting of 1,000 agents, which have opin- 
ion vectors of two dimensions. The normalization of 
Eq. © allows us to focus only on the temporal behav- 
ior of one of each agent's opinion without loss of in- 
formation. Fig. [5] exhibits this time behavior for the 
noises 77 = 0.05,0.2,0.35,0.45 during the first 100,000 
time steps. Each agent's first opinion is depicted by a 
color (gray tone) which corresponds to its value, and 
evolves beginning at the bottom. At low noise values 
stable clusters seem to form. The size of the clusters 
becomes smaller with decreasing 77. Nevertheless, these 
clusters are not stable, and the system reaches the con- 
sensus state after a finite time. For 77 = 0.05 and 77 = 0.2 
the size of clusters with the second dominant opinion is 
larger, indicating that at the end this opinion will con- 
trol the system. The larger the size of a cluster the longer 
it takes to break it. At values of 77 larger than 77 = 0.3, 
strong fluctuations control the system, and consensus be- 
gins to become unstable. For values around 77 = 0.3 one 
opinion still dominates, and clusters appear and disap- 
pear. At larger 77 the opinions have values around 0.5 for 
all agents which do not fluctuate much. 

It is interesting to calculate the number of time steps 
the system needs to reach its final state. In a system of 
two opinions we carried out various simulations with the 
same value of noise, 77, and a population size of N = 100. 
Each simulation begins with an initialization of randomly 
distributed opinions but a different random seed. Fig. [5] 
shows the distribution of times, needed to reach the con- 
sensus state for a system. The distribution decreases 
exponentially. The distribution becomes broader with 
decreasing values of 77, where 77 = should correspond to 
a flat distribution. 

As in the case of random interactions, a transition oc- 
curs from the consensus state to the one of a uniform 
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FIG. 5: (Color online) Graphical illustration of the temporal behavior of the system on a chain. The color (gray tone) 
corresponds to the value of the first of the two opinions of the system. The simulations run over 100, 000 time steps, drawing 
each 1, 000 iterations a new point on the vertical axis beginning at the bottom. The horizontal axis depicts the location of each 
agent on the chain, altogether consisting of 1,000 agents. The noise r\ is 0.05 in (a), 0.2 in (b), 0.35 in (c), and 0.45 in (d). 
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FIG. 6: Histogram of the time steps needed to reach the con- 
sensus state. Each curve corresponds to simulations with the 
same parameters: 100 agents, 2 opinions, random initializa- 
tion. For each value of the noise parameter we carried out 
200, 000 runs with different random seeds. 



distribution of the opinions. Fig. [7^ illustrates the vari- 



ation of (D) with r\ for different system sizes. As shown 
in Fig [7t>, by performing a finite-size scaling analysis 
through Eq. the collapse of all curves is obtained 
when we use the critical exponents v — 2, and [3 = 2. 



C. Complex networks 

In this section we compare the behavior of the opinion 
model if the agents interact with their k nearest neighbors 
on different networks topologies. We study two different 
kinds of scale-free networks; i.e. networks with a power 
law degree distribution k~ a . Those are the Barabasi- 
Albert network (BA) [l[ and the Apollonian network [l^] • 
The networks have considerable topological differences, 
that can be expressed in terms of their clustering coeffi- 
cient C. This coefficient is the average probability that 
the neighbors of a node are connected among them. The 
BA network has a clustering coefficient, C, which de- 
pends on the network size as iV _1 . It is independent 
on the degree of the nodes. In contrast, the Apollonian 
network has hierarchical structure with C depending on 
the degree of the node as a power law of the degree and 
its average value is high (C w 0.8) and independent of 
the network size N. Both types of scale-free networks, 
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FIG. 7: a: The fraction of agents with the same dominant 
opinion versus 77 are compared for different population sizes, 
N, in a system of two opinions, where only the nearest neigh- 
bors on a chain interact, b: Finite size scaling of the transi- 
tion. 



with and without hierarchical structure, have shown to 
be good models for rather different kinds of social inter- 
action networks, from social collaboration networks fl3j 
to networks of sexual contactsflij]. 

Further, we show that despite of the structural dif- 
ferences of these networks, the formation of consensus 
depends mainly on the noise and is independent on the 
specific topology of the scale-free network studied in the 
case of two opinions. The transition to consensus as a 
function of noise for the two scale-free networks, seems 
to belong to the same type of transition as in the case of 
annealed interactions. In contrast, we compare the be- 
havior of the model with a regular network with k = 6 
on a chain (in the previous section we had k = 2), adding 
interactions up to the third nearest neighbors. The tran- 
sition from consensus to a uniform distribution on the 
regular network differs from the transition of complex 
networks and annealed case and presents similar behav- 
ior as the one reported in previous sections for a chain. 

In Fig. [H^ we show (D) vs. 77 for the model on the BA 
(triangles) and Apollonian networks (solid line), com- 
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FIG. 8: (Color online) a:Influence of the topology of networks 
on the transition to consensus (D — 1) as a function of noise 
(77). The transition on two different scale- free networks, the 
Apollonian (solid line) and the Barabasi- Albert network (BA, 
triangles) is similar to the one observed for annealed interac- 
tions (AI, plus signs), and differs from the transition on a 
regular lattice (circles). In the three insets we plot the value 
of the dominant opinion, a 0max vs. time, b: Comparison 
of the behavior of a 0rnax (t) on the four networks: Apollonian 
(solid line), Barabasi- Albert (BA, dotted line), Annealed (AI, 
dashed-line) and Regular (long dashed line) for a fixed noise 
(77 = 0.2). One observes that for this noise, which is below the 
critical noise, in the regular network the emergence of consen- 
sus takes longer than in scale-free and annealed interactions, 
which have similar behavior (three upper curves), c: Near 
below the transition, for r\ = 0.4, we compare the response 
of the regular and the Apollonian network. It is observed 
that for the former there is an intermittency among consen- 
sus D = 1 and a 0max = 0.79 and not consensus D — 0.5 and 
a 0mam = 0.5. This behavior is not observed in the complex 
networks, d: Above the transition (77 = 0.6) the consensus is 
broken and the dynamics of the opinion a 0rnax vs. time be- 
haves similar in regular and complex networks. All simulation 
runs are with systems of 124 agents. 



pared to the result of annealed interactions (plus signs) 
and the regular network (circles). The results of the fig- 
ure represent the average over 20 realizations on systems 
of N = 124 agents and 2 opinions. Near the transition, 
the fluctuations on the regular lattice strongly increase, 
as opposed to annealed interaction and to BA and Apol- 
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Ionian networks. This is because the system presents an 
intermittency near the transition point (77 rj 0.4). We 
observe this intermittency of the dynamics in Fig. [8fc, 
comparing the value of the dominant opinion a 0max vs. 
time, for the Apollonian and the regular network with 
77 = 0.4. Above the critical noise, there is no consen- 
sus and the fraction of agents with dominant opinion is 
(D) Ri 1/2. At these values of 77, the response of the sys- 
tem is similar for scale-free and regular networks, as is 
shown in Fig. [5U with 77 = 0.6. 

Above the critical noise there is no way for the agents 
to achieve global coordination. In this situation, the dy- 
namics is dominated by local interactions, thus the topol- 
ogy of the system has little effect on this regime. Be- 
low the critical noise, global coordination becomes pos- 
sible. However, the low dimensionality of the regular 
lattice leads to the intermittent behavior observed in the 
panel [8b. 

IV. CONCLUSION 

Starting from a model based on interactions with a 
term of second order, we analyzed its behavior for dif- 
ferent topologies: random, regular and complex ones. 
Depending on the control parameter, the noise 77, two 
different absorbing states control the system. Its behav- 
ior changes from consensus to a uniform distribution of 
opinions. Despite the rather simple approach to take 
into account such simple interactions, a rich variety of 
results can be reported depending on the dimension of 



the opinion vector. The results show that an opinion is 
kept (for systems with more than two opinions), and the 
parameters need to be adjusted crucially to change the 
state (hysteresis). This occurs at different dimensions 
O of the opinion vector, depending on the topology of 
interactions. 

The response of the system to approach consensus has 
the origin in the model dynamics as opposed to the par- 
ticular features of the network. An important charac- 
teristic of the transition to consensus is the dimension 
associated to the space of agent interactions. The dy- 
namical response of the opinion model for both scale-free 
networks is similar to the response observed for annealed 
interactions and each of these cases represents long range 
interactions. In contrast, differences are reported with a 
regular lattice, which has spatial dimension one, associ- 
ated to nearest neighbors interactions. 

As it was previously observed for the Sznajd model of 
opinion formation, for the general model that we present 
here, the response of the system in terms of opinion for- 
mation is qualitatively the same for a deterministic scale 
free network as well as for a random scale free network. 
This implies a clear advantage for an analytical treatment 
in a similar way as was done for the Sznajd model [15j. 
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